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HIGH-SPEED INVERSE DISCRETE COSINE TRANSFORMATION 
METHOD AND APPARATUS 

BACKGROUND OF THE INVENTION 

[01] This application claims the priority of Korean Patent Application No. 
2002-72384, filed on November 20, 2002, in the Korean Intellectual Property 
Office, the disclosure of which is incorporated herein in its entirety by 
reference. 

1 . Field of the Invention 

[02] The present invention relates to a high-speed inverse discrete cosine 
transformation method and apparatus. 

2. Description of the Related Art 

[03] Compression of digital data, such as image signals, is one of the most 
important techniques in an environment that supports multimedia applications. 
Since image signals generally consist of a considerable amount of data, 
effectively transmitting, storing, and processing image signals has been 
limited. In order to overcome such restrictions, numerous compression stream 
grammars and decoding techniques have been proposed through various 
international standards, such as MPEG-2, MPEG-4, H.263, and H.26L. 

[04] There are two different types of compression techniques, i.e., a lossless 
compression technique and a loss compression technique. When adopting the 
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lossless compression technique, data, such as characters, figures, or other 
ordinary data, can be compressed at an average compression rate of 2:1, and 
the compressed data can be flawlessly restored. On the other hand, when 
adopting a loss-compression technique, i.e., when compressing image data, 
voice data, or acoustic data, minor data loss that is imperceptible to a person is 
allowed and a compression rate of 10:1 can be achieved. One of the most 
common loss compression techniques is conversion encoding. In conversion 
encoding, data is arranged in a predetermined manner having high spatial 
correlation with one another and subjected to orthogonal conversion. During 
orthogonal conversion, data is divided into a variety of frequency components 
ranging from low-frequency components to high-frequency components, and 
then each of the frequency components is quantized. By doing so, the 
correlation among the frequency components almost disappears, and signal 
energy is concentrated at a low-frequency range. Among the low-frequency 
components resulting from orthogonal conversion, the components on which 
more energy is concentrated, i.e., the ones having a higher dispersion value, 
are more accurately represented through the use of additional bits. A low- 
frequency component having a dispersion value four times greater than the 
dispersion values of other components (i.e., a low-frequency component 
having an amplitude two times greater than the amplitude of other 
components) is assigned an additional bit. Finally, all the frequency 
components are expected to have the same quantization error characteristics. 
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[05] Among the various types of orthogonal conversion, Karhunen-Loeve 
Transformation is considered one of the most effective compression 
techniques because image signals subjected to this transformation technique 
have superior energy concentration characteristics. However, Karhunen- 
Loeve Transformation requires different conversion functions for different 
images, which imposes serious restrictions on the implementation of 
Karhunen-Loeve Transformation. As an alternative to Karhunen-Loeve 
Transformation, which is difficult to apply, discrete cosine transformation 
(DCT) has been suggested. Since DCT exhibits almost the same performance 
results as Karhunen-Loeve Transformation, can be practically applied, and is 
also realizable, it is considered one of the core technologies in a variety of 
international standards. In the DCT technique, 8x8 pixels are grouped into 
one block and then each block is discrete-cosine- transformed. As the size of 
blocks increases, the efficiency of data compression becomes higher, but it 
becomes more difficult to perform DCT on each block. Through a number of 
experiments, an 8 x 8 block has been determined as the DCT unit which can 
meet both the requirements of efficient data compression and easy 
implementation. 

[06] Conventional data compression techniques have used discrete cosine 
transformation to eliminate spatial redundancies that are obtained when 
compressing images. Motion estimation (ME) and motion compensation 
(MC) have been used to eliminate temporal redundancies. 
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[07] FIG. 1 is a diagram illustrating conventional two dimensional inverse 
discrete cosine transformation (2D-IDCT). High-speed IDCT is achieved by 
reducing the amount of computations necessary for IDCT. As shown in FIG. 
1, 2D IDCT is characterized as performing one-dimensional (ID) IDCT on 
each row and each column of a block. In particular, the 2D IDCT shown in 
FIG. 1 can be characterized by ID IDCT performed along two different 
directions, displacement, and a row column algorithm. In other words, in 2D 
IDCT, 1 D IDCT is performed on a block N times in a row direction, the block 
is displaced, and then IDCT is performed on the displaced block N times in a 
column direction. 

[08] FIG. 2 is a table for comparing the number of computations for 
different conventional IDCT algorithms. There are four different conventional 
IDCT algorithms used in the 2D IDCT of FIG. 1, i.e., Chen's algorithm based 
upon sparse matrix factorizations, Wang's algorithm, Lee's algorithm based 
upon discrete Fourier transformation, and the AAN algorithm. FIG. 2 shows 
the number of multiplications and additions required for each of the 
conventional IDCT algorithms. 

[09] FIG. 3 is a diagram illustrating Wang's algorithm, which is one of the 
high-speed IDCT algorithms. Referring to FIG. 3, signals X(0) through X(7) 
are output in response to input signals x(0) through x(7), respectively. In the 
process of outputting the signals X(0) through X(7), the calculation result 
carried out on a contact point must be passed on to a following contact point. 
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A contact point can be obtained by multiplying the previous contact point by a 
constant C 0 through C 7 placed therebetween. For example, C 0 = x(l) + x(7). 
Wang's algorithm needs 16 adders and 26 multipliers. Most high-speed 
algorithms are more or less similar to Wang's algorithm mentioned above. 

[10] Conventional high-speed IDCT algorithms can generally reduce the 
complexity of computations necessary for IDCT. In the process of restoring 
compressed data, however, conventional high-speed IDCT algorithms require 
a considerable number of computations. In the current mobile environment, 
which is capable of providing a variety of multimedia services, decoders, i.e., 
mobile communications devices such as mobile phones or personal digital 
assistants (PDAs), are restricted in terms of size and power consumption, 
while encoders, i.e., multimedia service providers' server systems, are 
relatively free from those restrictions. Therefore, there is a need to reduce the 
amount of computations necessary for performing IDCT in decoders. 

SUMMARY OF THE INVENTION 

[11] The present invention provides a high-speed inverse discrete cosine 
transformation (IDCT) method and apparatus, which are capable of 
considerably reducing the number of computations during IDCT by 
performing two-dimensional (2D) IDCT on a discrete cosine transformation 
(DCT) matrix on an element-by-element basis or on a matrix-by-matrix basis, 
depending on the number of elements with an valid value. 
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[12] According to an aspect of the present invention, there is provided a 
high-speed inverse discrete cosine transformation (IDCT) method, which 
involves (a) searching all elements of a discrete cosine transformation (DCT) 
matrix for elements having a value other than 0, in a predetermined order, 
when a total number of elements having a value other than 0 is not greater 
than a predetermined critical value; (b) performing two-dimensional (2D) 
IDCT on the elements having a value other than 0 searched in (a); and (c) 
performing 2D IDCT on the DCT matrix when the total number of elements 
having a value other than 0 is greater than the predetermined critical value. 

[13] According to another aspect of the present invention, there is provided 
a high-speed IDCT apparatus, including an element searching unit, an 
element-wise 2D IDCT unit, and a matrix-wise 2D IDCT unit. The element 
searching unit searches all elements of a discrete cosine transformation (DCT) 
matrix for elements having a value other than 0 in a predetermined order, 
when a total number of elements having a value other than 0 is not greater 
than a predetermined critical value. The element-wise 2D IDCT unit performs 
2D IDCT on the elements having a value other than 0 searched by the element 
searching unit. The matrix-wise 2D IDCT unit performs 2D IDCT on the 
DCT matrix when the total number of elements having a value other than 0 is 
greater than the predetermined critical value. 
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BRIEF DESCRIPTION OF THE DRAWINGS 
[14] The above and other features and advantages of the present invention 
will become more apparent by describing in detail exemplary embodiments 
thereof with reference to the attached drawings in which: 

[15] FIG. 1 is a diagram illustrating a conventional two-dimensional inverse 
discrete cosine transformation; 

[16] FIG. 2 is a table showing the number of computations for different 
conventional inverse discrete cosine transformation (IDCT) algorithms; 

[17] FIG. 3 is a diagram illustrating Wang's algorithm, which is one of the 
conventional high-speed IDCT algorithms; 

[18] FIG. 4 is a block diagram of a high-speed IDCT apparatus according to 
a preferred embodiment of the present invention; 

[19] FIG. 5 is a diagram illustrating an example of a typical 8x8 discrete 
cosine transformation (DCT) matrix; 

[20] FIG. 6 is a block diagram of an element-wise 2D IDCT unit of FIG. 4; 

[21] FIG. 7 is a flowchart of a high-speed IDCT method according to a 
preferred embodiment of the present invention; and 

[22] FIG. 8 is a table for comparing the number of computations used in 
high-speed IDCT according to the present invention to the number of 
computations used in conventional high-speed IDCT. 
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DETAILED DESCRIPTION OF THE INVENTION 

[23] Hereinafter, the present invention will be described in greater detail 
with reference to the accompanying drawings in which preferred embodiments 
of the invention are shown. 

[24] FIG. 4 is a block diagram of a high-speed inverse discrete cosine 
transformation (IDCT) apparatus according to a preferred embodiment of the 
present invention. Referring to FIG. 4, the high-speed IDCT apparatus 
includes an effective element number calculation unit 40, an element 
searching unit 41, an element- wise two-dimensional (2D) IDCT unit 42, and a 
matrix-wise 2D IDCT unit 43. 

[25] A process for compressing image data is as follows. First, an image 
signal is discrete-cosine-transformed so that the image signal is divided into 
several frequency ranges. Here, the energy of image data is generally 
concentrated in low- frequency ranges. Therefore, by quantizing image data, it 
is possible to compress the image data using a reduced number of bits. 
Thereafter, quantization is performed on the image data so that the image data 
is divided into identical sized quanta. Then, quanta represented by a value 
smaller than a predetermined number are given a value of 0 to replace their 
respective original data values, so that the size of the entire data can be 
reduced. At this moment, data loss may occur due to the assignment of 0 
values. Thereafter, run-length encoding is performed, in which a repetition of 
characters is replaced with the number of same characters and a single 
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character. As the number of characters constituting each run increases and the 
number of run occurrences increases, compression efficiency increases. 
Thereafter, Huffman encoding is performed, in which integer sequences 
obtained by zigzag scanning are converted into binary values. By doing so, an 
8x8 matrix is compressed into several combinations of Os and Is. In order to 
decode such compressed image data, the above-described compression process 
must be inversely performed. 

[26] As described above, in the case of compressing image data according 
to the process of performing DCT and quantization on the image data, the 
compressed data is mostly concentrated in low-frequency ranges, while almost 
nothing is left in high-frequency ranges. As the quantization scale becomes 
larger, the number of elements having a value of 0 increases while the number 
of elements having a value other than 0 decreases. In contrast, as the 
quantization scale becomes smaller, the number of elements having a value 
other than 0 increases. In a scenario with a large quantization scale, the 
number of elements in a DCT matrix having a value other than 0 is small. 
Using a conventional high-speed IDCT algorithm, only these elements are 
inversely discrete-cosine- transformed using a conventional high-speed IDCT 
algorithm without the need to inversely discrete-cosine-transform all elements 
of the DCT matrix. On the other hand, in a scenario with a small quantization 
scale, many elements having a value other than 0 exist, and it is effective to 
use the conventional high-speed IDCT algorithm. 
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[27] In a case where the number of elements in the DCT matrix having a 
value other than 0 is not greater than a predetermined critical value, the 
element searching unit 41 searches for the elements having a value other than 
0 in a predetermined order. The element-wise 2D IDCT unit 42 performs 2D 
IDCT on the searched elements. As described above, the number of elements 
having a value other than 0 is compared with the predetermined critical value 
to determine whether IDCT is performed on an element-by-element basis or 
on a matrix-by-matrix basis. The predetermined critical value represents a 
maximum number of elements having a value other than 0, at which the 
number of computations for element-wise IDCT is less than the number of 
computations for matrix-wise IDCT. A number of experiments has 
determined that with a quantization scale not greater than 10, the highest 
compression efficiency is achieved when the critical value is set to 15, and 
with a quantization scale greater than 1 0, the highest compression efficiency is 
achieved when the critical value is set to 10. 

[28] As described above, the process of decoding a compressed file is just 
the opposite of the process of encoding a file into the compressed file. 
Therefore, the decoding process is carried out by sequentially performing 
Huffman decoding, run-length decoding, inverse quantization, and IDCT on 
the compressed file. The number of elements in the DCT matrix having a 
value other than 0 can be determined in advance during the process of run- 
length decoding, which is performed prior to IDCT. Since run-length 
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encoding replaces a series of Os with a single 0 and the length of the series of 
Os, the length of the series of Os can be figured out in run-length decoding. 
Through run-length decoding it is also possible to identify the number of 
elements having a value other than 0. 

[29] In other words, the effective element number calculation unit 40 counts 
only elements having a value other than 0 in the process of run-length 
decoding. For example, suppose that a quantization scale is not greater than 
10 and a critical value is set to 15 in the process of encoding, and the number 
of elements having a value other than 0 is not greater than 15, it is effective to 
perform IDCT on an element-by-element basis rather than on a matrix-by- 
matrix basis. 

[30] In a case where the number of elements in the DCT matrix having a 
value other than 0 is greater than the predetermined critical value, the matrix- 
wise 2D IDCT unit 43 performs 2D IDCT on the DCT matrix. In this case 
where a quantization scale is not greater than 10 and a critical value is set to 
15 in the process of encoding, and the number of elements having a value 
other than 0 is greater than 15, it is effective to perform IDCT on a matrix-by- 
matrix basis rather than on an element-by-element basis. The matrix-wise 2D 
IDCT unit 43 performs such as 2D IDCT on the DCT matrix using a 
conventional high-speed IDCT algorithm such as Wang's algorithm, Chen's 
algorithm, Lee's algorithm, or the AAN algorithm. 
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[31] FIG. 5 is a diagram illustrating a typical example of an 8x8 DCT 
matrix. In FIG. 5, an element representing more concentrated image data is 
darker. Element Al, which belongs to the lowest frequency range, is the 
darkest one in the 8x8 DCT matrix. In the 8x8 DCT matrix, the value of 
element Al is called a direct current (DC) value, and values of the other 63 
elements, A2 through A64, are called alternating current (AC) values. The 
DC value and the values of elements neighboring element Al are very 
important in determining the brightness of the entire image. 

[32] As for the predetermined order used by the element searching unit 41 
when searching for the elements having a value other than 0, the element 
searching unit 41 searches the DCT matrix along a horizontal direction or a 
vertical direction. Since most of the elements having a value other than 0 are 
concentrated near the upper left corner of the DCT matrix, it is preferable for 
the element searching unit 41 to search the elements of the DCT matrix in a 
zigzag manner, i.e., in the order of Al, A2, A3, ... , A64, thereby reducing the 
number of calculations. 

[33] FIG. 6 is a block diagram of the element-wise 2D IDCT unit 42 of 
FIG. 4. Referring to FIG. 6, the 2D IDCT unit 42 includes a partial value 
calculator 61 and a complete value calculator 62. 

[34] In general, 2D IDCT is performed using Equation (1) below. 
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ifi±0andj*0, C(iJ)= — 
ifi=0orj=0, C(i 9 j)=jj 

[35] In Equation (1), 7(7, j) represents the value of an element located at 
(i+1, j+1) of a DCT matrix T, and V(x, y) represents the value of an element 
located at (x+1, y+1) of a matrix V, which represents the restored matrix 
obtained through IDCT on the DCT matrix T. If the DCT matrix T is an 8x8 
matrix, ij 9 jt, and y each have a value between 0 and 7, and N=8. 

[36] A superimposition principle is adopted to selectively process the 
elements in the DCT matrix having a value other than 0. According to the 
superimposition principle, each DCT coefficient block or every predetermined 
number of DCT coefficient blocks is inversely discrete-cosine-transformed, 
and then all IDCT results are summed up. This process achieves the same 
results as those produced after performing IDCT on all the DCT coefficient 
blocks of the DCT matrix at the same time. In the present invention, DCT 
coefficient block values are searched one by one, DCT coefficient blocks 
having valid values are inversely discrete-cosine-transformed, and IDCT 
results are summed up, thus obtaining a restored matrix (restoring an original 
version of the DCT matrix). In order to perform 2D IDCT on an element-by- 
element basis, the predetermined calculation process shown in Equation (2) 
must be performed on elements in the DCT matrix having a value other than 0. 
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IDCT(T)= IDCT(T(0,0))+ IDCT(T(0,l)) + IDCT(T(0,2))+--+ IDCT(T(OJ)) 
+ /DC7X 7X1,0)) + IDCT(T(U))+ IDCT(T{\2))+ IDCT(T(IJ)) 
+ /Z>C7X 7X2,0))+ 7DCr(r(2,l))+ /DCr(7 T (2,2))+ - -+ IDCT(T(2J)) 
+ •••+ 

+ IDCT(T(7,0))+ 7DCr(r(7,l))+ IDCT(T(JZ))+ • + IDCT(T(7 J)) 

[37] In Equation (2), IDCT(T) represents an 8x8 matrix V restored from the 
8x8 DCT matrix T. IDCT(T(0, 0)) represents an 8x8 matrix, which is 
obtained by substituting the variables in Equation (1) with the coordinates (i=0 
and j=0) and value T(0, 0) of the element in the first row and the first column 
of the DCT matrix T, and the coordinates (x and y are a value between 0 and 
7) of each element in the restored 8x8 matrix V. IDCT(T(0, 0)) accounts for 
part of the restored 8x8 matrix V. In other words, according to the 
superimposition principle, the restored matrix V is obtained by summing up all 
matrices generated for the elements having a value other than 0. In the present 
invention, each of the matrices generated for the elements having a value other 
than 0 is stored in memory in a table format, and values located at memory 
addresses and corresponding to all the tables stored in the memory are 
summed up, thereby obtaining the restored matrix V. 

[38] The above-described algorithm can be applied to the case where the 
number of elements having a value other than 0 is not greater than the critical 
value, given that the critical value, which is used to determine whether to 
perform IDCT on an element-by-element basis or on a matrix-by-matrix basis, 
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can be 6, 10, or 1 5, depending on the quantization scale. In the case where the 
number of elements having a value other than 0 is greater than the critical 
value, a conventional high-speed algorithm can be used. Various image 
compression algorithms, such as MPEG-2, MPEG-4, and H.261, generally 
produce no more than 10 elements having a value other than 0. 

[39] For example, when the number of elements having a value other than 0 
is 10, element-wise 2D IDCT is represented by Equation (3) below. 

IDCT(T)= IDCT(T(0,0))+ IDCT( 7(0,1)) + IDCT( 7(0,2)) + /DC7(7(0,3)) — (3) 
= IDCT( 7(1,0)) + IDCT( 2TU))+ IDCT( 7(1,2)) 
= IDCT( 7(2,0))+ /£>C7( 7(2,1)) 
= /£>C7(7(3,0)) 

[40] The partial value calculator 61 obtains IDCT(T(i,j)), which constitutes 
part of the restored matrix V, by substituting variables in Equation (1) with the 
coordinates of an element having a value other than 0, which is searched 
for by the element searching unit 41 of FIG. 4, and the coordinates (x, y) of 
each element in the restored matrix V. The complete value calculator 62 takes 
advantage of the superimposition principle and obtains a complete value V(x, 
y) corresponding to each element of the restored matrix V. The value is 
obtained by summing up IDCT(T(iJ)), which is generated by the partial value 
calculator 61 for each of the elements that has a value other than 0 and has 
been searched for by the element searching unit 41. 

[41] FIG. 7 is a flowchart of a high-speed IDCT method according to a 
preferred embodiment of the present invention. Referring to FIG. 7, in step 
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70, the number of elements in a DCT matrix having a value other than 0 is 
counted during run-length decoding, which is part of a process of the decoding 
process for a predetermined compressed file and is performed before IDCT. 

[42] If the number of elements having a value other than 0 is not greater 
than a predetermined critical value in step 71, the elements having a value 
other than 0 are searched for one by one in a predetermined order in step 72. 
Here, the number of elements having a value other than 0 is discovered in the 
process of run-length decoding, which is performed before IDCT. In addition, 
the critical value is set to the maximum number of elements having a value 
other than 0, at which the number of computations for element-wise IDCT is 
less than the number of computations for matrix-wise IDCT. Furthermore, the 
search for elements having a value other than 0 include a zigzag search across 
all elements of the DCT matrix, starting with the element in the first column 
and first row of the DCT matrix. 

[43] Thereafter, the elements having a value other than 0 are inversely 
discrete-cosine- transformed in step 73. In other words, a partial value for each 
element of a restored matrix is obtained by substituting the variables in an 
IDCT formula with the value and coordinates of each element having a value 
other than 0 and the coordinates of each element in the restored matrix. 
Thereafter, according to the superimposition principle, partial values obtained 
for all the elements of the restored matrix are summed up, thereby obtaining a 
complete value for each element of the restored matrix. 
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[44] If the number of elements having a value other than 0 in the DCT 
matrix is greater than the critical value in step 71, the DCT matrix is subjected 
to 2D IDCT in step 74. In other words, the DCT matrix is two-dimensionally 
inverse-discrete-cosine-transformed using a conventional high-speed IDCT 
algorithm, such as Wang's algorithm, Chen's algorithm, Lee's algorithm, or 
the AAN algorithm. 

[45] FIG. 8 is a table for comparing the number of computations necessary 
for high-speed IDCT according to the present invention to the number of 
computations necessary for conventional high-speed IDCT. FIG. 8 shows 
experimental results obtained when a quantization scale is large. As described 
above, in a case of a small quantization scale, use of a conventional 2D IDCT 
algorithm is required. 

[46] The above-described embodiments of the present invention can be 
realized as a computer program that can be recorded on a computer-readable 
recording medium and can be executed in a digital computer. 

[47] The computer-readable recording medium includes a magnetic storage 
medium, such as ROM, a floppy disk, or a hard disk; an optical recording 
medium, such as a CD-ROM or a DVD; and a carrier wave, such as data 
transmission through the Internet. 

[48] According to the present invention, it is possible to minimize the 
number of computations by eliminating unnecessary computations for 
elements having a value of 0, which account for the majority of DCT matrix 
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elements. In addition, the present invention provides an optimum IDCT 
algorithm for different quantization scales. For example, when the data 
compression rate is low and the number of elements having a value other than 
0 is small, a conventional high-speed IDCT algorithm is adopted. Even 
though the element-wise 2D IDCT algorithm of the present invention and a 
conventional IDCT algorithm are both used, the number of computations for 
IDCT are still considerably reduced, because nearly 80% of target image 
signals are subjected to the element- wise 2D IDCT algorithm of the present 
invention. The percentage of image data that can be processed using the 
element-wise 2D IDCT algorithm of the present invention, rather than the 
conventional high-speed IDCT algorithm, varies on a case-by-case basis. 
Furthermore, according to the present invention, it is possible to design a 
stable video decoder having enhanced performance or a compact-sized mobile 
video decoder having reduced power consumption by dramatically reducing 
the number of computations performed by an IDCT module, which amounts to 
nearly 25 - 30 % of the total number of computations performed in a video 
decoder. 

[49] In the prior art, computations in high-speed DCT algorithms are 
carried out with reference to the end of block (EOB) and different high-speed 
IDCT algorithms are required for different scanning methods. However, the 
present invention can be applied irrespective of the type of scanning, whether 
the type is zigzag scanning, horizontal scanning, or vertical scanning. In 
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addition, the present invention maximizes its use of valid numbers in 
computations. In other words, since in the present invention, only one round 
of IDCT is carried out, the peak signal-to-noise ratio is higher in the present 
invention than in the prior art. 

[50) While the present invention has been particularly shown and described 
with reference to exemplary embodiments thereof, it will be understood by 
those of ordinary skill in the art that various changes in form and details may 
be made therein without departing from the spirit and scope of the present 
invention as defined by the following claims. 
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